On paired decoupled quasi-linearization methods for solving nonlinear systems of differential equations that model boundary layer fluid flow problems.

Doctoral Degree. University of KwaZulu-Natal, Pietermaritzburg.Two numerical methods, namely the spectral quasilinearization method (SQLM) and the spectral local linearization method (SLLM), have been found to be highly efficient methods for solving boundary layer flow problems that are modeled using systems of differential equations. Conclusions have been drawn that the SLLM gives highly accurate results but requires more iterations than the SQLM to converge to a consistent solution. This leads to the problem of figuring out how to improve on the rate of convergence of the SLLM while maintaining its high accuracy. The objective of this thesis is to introduce a method that makes use of quasilinearization in pairs of equations to decouple large systems of differential equations. This numerical method, hereinafter called the paired quasilinearization method (PQLM) seeks to break down a large coupled nonlinear system of differential equations into smaller linearized pairs of equations. We describe the numerical algorithm for general systems of both ordinary and partial differential equations. We also describe the implementation of spectral methods to our respective numerical algorithms. We use MATHEMATICA to carry out the numerical analysis of the PQLM throughout the thesis and MATLAB for investigating the influence of various parameters on the flow profiles in Chapters 4, 5 and 6. We begin the thesis by defining the various terminologies, processes and methods that are applied throughout the course of the study. We apply the proposed paired methods to systems of ordinary and partial differential equations that model boundary layer flow problems. A comparative study is carried out on the different possible combinations made for each example in order to determine the most suitable pairing needed to generate the most accurate solutions. We test convergence speed using the infinity norm of solution error. We also test their accuracies by using the infinity norm of the residual errors. We also compare our method to the SLLM to investigate if we have successfully improved the convergence of the SLLM while maintaining its accuracy level. Influence of various parameters on fluid flow is also investigated and the results obtained show that the paired quasilinearization method (PQLM) is an efficient and accurate method for solving boundary layer flow problems. It is also observed that a small number of grid-points are needed to produce convergent numerical solutions using the PQLM when compared to methods like the finite difference method, finite element method and finite volume method, among others. The key finding is that the PQLM improves on the rate of convergence of the SLLM in general. It is also discovered that the pairings with the most nonlinearities give the best rate of convergence and accuracy

On decoupled quasi-linearization methods for solving systems of nonlinear boundary value problems.

M. Sc. University of KwaZulu-Natal, Pietermaritzburg 2014.In this dissertation, a comparative study is carried out on three spectral based numerical methods which are the spectral quasilinearization method (SQLM), the spectral relaxation method (SRM) and the spectral local linearization method (SLLM). The study is carried out by applying the numerical methods on systems of differential equations modeling uid ow problems. Residual error analysis is used in determining the speed of convergence, convergence rate and accuracy of the methods. In Chapter 1, all the terminologies and methods that are applied throughout the course of the study are introduced. In Chapter 2, the SRM, SLLM and SQLM are applied on an unsteady free convective heat and mass transfer on a stretching surface in a porous medium with suction/injection. In Chapter 3, the SRM, SLLM and SQLM are applied on an unsteady boundary layer ow due to a stretching surface in a rotating uid. In Chapter 4, the SRM, SLLM and SQLM are used to solve an unsteady three-dimensional MHD boundary layer ow and heat transfer over an impulsively stretching plate. The purpose of this study is to assess the performance of the spectral based numerical methods when solving systems of differential equations. The performance of the methods are measured in terms of computational efficiency (in terms of time taken to generate solutions), accuracy and rate of convergence. The ease of development and implementation of the associated numerical algorithms are also considered